Aharonov-Bohm effect for a valley-polarized current in graphene 
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This is a numerical study of the conductance of an Aharonov-Bohm interferometer in a tight- 
binding model of graphene. Two single-mode ballistic point contacts with zigzag edges are connected 
by two arms of a hexagonal ring enclosing a magnetic flux $. The point contacts function as valley 
filters, transmitting electrons from one valley of the band structure and reflecting electrons from 
the other valley. We find, in the wider rings, that the magnetoconductance oscillations with the 
fundamental periodicity A$ = h/e are suppressed when the two valley filters have opposite polarity, 
while the second and higher harmonics are unaffected or enhanced. This frequency doubling is 
interpreted in terms of a larger probability of intervalley scattering for electrons that travel several 
times around the ring. In the narrowest rings the current is blocked for any polarity of the valley 
filters, with small, nearly sinusoidal magnetoconductance oscillations. Qualitatively similar results 
are obtained if the hexagonal ring is replaced by a ring with an irregular boundary. 

PACS numbers: 73.23.-b, 73. 23. Ad, 73.63.Rt, 85.35.Ds 



I. INTRODUCTION 

The isolation of single two-dimensional layers of carbon 
(graphene) [l| has led experimental and theoretical physi- 
cists to reexamine classic effects from mesoscopic physics 
0. The search is for novel features that arise from the 
unusual conical band structure of a carbon monolayer Q . 
One effect which has so far received little attention is the 
periodic magnetoconductance oscillation in a ring known 
as the Aharonov-Bohm effect [1, Q . Preliminary experi- 
ments on the magnetoconductance of graphene rings have 
been reported by several groups [1, 01 ■ Theoretically, the 
energy spectrum of a closed graphene ring was studied 
as a function of the enclosed magnetic flux $ by Recher 
et al. [|. 

Here we study by computer simulation the electron 
transport through an open graphene ring, contacted to 
electron reservoirs by ballistic point contacts. Our work 
builds on an earlier finding [qJ that a single-mode point 
contact with zigzag edges operates as a valley filter. De- 
pending on whether the Fermi level in the point contact 
lies in the conduction or valence band, the transmitted 
electrons occupy states in one or the other valley of the 
band structure. It was also shown in Ref. that two 
adjacent valley filters function as a highly effective val- 
ley valve, passing or blocking the current depending on 
whether the two filters have the same or opposite polar- 
ity. 

We find that the magnetoconductance oscillations in a 
multimode graphene ring with valley-filtering point con- 
tacts show the expected A$ = h/e periodicity if the two 
filters have the same polarity. For opposite polarity, how- 
ever, a period doubling appears: The lowest harmonic is 
suppressed while the second and higher harmonics are 
unaffected or enhanced. We attribute the period dou- 
bling to intervalley scattering, which is more effective for 



electrons that have travelled more than once along the 
ring. 

In few-mode rings the period doubling does not hap- 
pen, instead the conductance is strongly suppressed with 
small, nearly sinusoidal, magnetoconductance oscilla- 
tions. They appear in their clearest form in a hexago- 
nal ring, but we find qualitatively the same behavior in 
a more generic ring with irregular boundaries. 



II. AHARONOV-BOHM INTERFEROMETER 
IN GRAPHENE 

The analysis starts from the tight-binding model of 
graphene, with Hamiltonian 

H = Y,n^^){J\+Y.n^)(^\■ (1) 

The system is coupled to the vector potential A through 
the hopping matrix element 

with T = 3 eV the hopping energy and = h/e the flux 
quantum. The orbitals \i) and \j) are nearest neighbors 
on a honeycomb lattice (with lattice points Ri), other- 
wise Tij = 0. The energy-independent velocity v near the 
Dirac point equals v = ^y/Sra/h « 10^ m/s, with the 
lattice constant a = 0.246 nm. 

The electrostatic potential Vi — V{xi) varies only 
along the axis connecting the input and output point 
contacts (see Fig. [T]). Namely, the potential equals U at 
the first constriction (0 < a; < Z, where / is the constric- 
tion length), U' at the second constriction {l + Lc < x < 
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FIG. 1: Schematic diagram of a hexagonal graphene ring at- 
tached to graphene leads (top panel) and corresponding po- 
tential profile (bottom panel). The constriction with zigzag 
edges between each lead and the ring allows one to control 
input and output valley polarizations by varying the electro- 
static potentials U and U' . 



at i? = 35-\/3a and vary the inner radius. We take ring 
widths W/^/^a = 5, 10, 15, corresponding to = 1, 3, 5. 
(Even values of N are not accessible because of the val- 
ley degeneracy of the second and higher modes.) The 
external leads have a fixed width Woo/^f^a = 70, corre- 
sponding to 29 propagating modes. 

We take the vector potential A = {A^-, 0, 0) with 
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< X < 2l + Lc- 
othcrwise. 
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This corresponds to a uniform perpendicular magnetic 
field B in the area containing the ring, the two point 
contacts, and the widening region connecting the point 
contacts to the external leads (see grey rectangle in Fig. 
[TJ top panel). For technical reasons, the magnetic field 
is set to zero in the external leads. The resulting step in 
the field strength at the entrance and exit of the interfer- 
ometer will reduce its conductance somewhat, but this 
is not expected to change the qualitative features of the 
magnetoconductance oscillations that we are interested 
in. 

We calculate the transmission matrix t numerically and 
then obtain the conductance from the Landauer formula 
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(The factor of two accounts for the spin degeneracy.) 



III. RESULTS 



21 + Lc: with Lc — (4i? — 'w)/\/3, where R is the ring 
radius and w is the constriction width), and zero every- 
where else. By varying U and U' at a fixed Fermi energy 
Hoo in the external leads, we can vary the Fermi energies 
A* = Moo — U and /i' — ^oo — U' in the two constrictions. 
We took Hoc = t/3 to work with heavily doped graphene 
leads, while remaining at sufficiently small Fermi energy 
that the linearity of the dispersion relation holds reason- 
ably well. The constriction parameters w = 10\/3a and 
I — 16a are chosen to provide valley polarizations above 
90% i. 

We denote the number of transmitted modes through 
the first constriction by n and through the second con- 
striction by n'. A positive number indicates that the 
Fermi level lies in the conduction band, while a nega- 
tive number indicates that the Fermi level lies in the va- 
lence band. For example, as shown in Ref. the case 
-3A/2 < /i' < < ^ < 3A/2 (with A = ttHv/w) cor- 
responds to n = 1, n' — —1 (valley filters of opposite 
polarity), while the case < < 3 A/2 corresponds 
to n = n' = 1 (valley filters of the same polarity) . 

The AB interferometer is modeled by a hexagon with 
a hexagonal hole in the center and zigzag edges along 
the entire perimeter. (We will consider a more generic 
shape in Sec. IIVI ) To vary the number of modes TV that 
can propagate along the ring, we keep the radius fixed 



To analyze the operation of the hexagonal graphene 
ring as a valley valve 0], we first show in Fig. 2 the 
conductance in zero magnetic field. We take /i = 0.05 r 
to keep n = 1 at the first constriction (the polarizer), 
and vary the Fermi energy /i' at the second constriction 
(the analyzer). For negative /i' {n' — —1), when the 
two constrictions transmit the opposite valley polariza- 
tion, we obtain G < 10"'^ e^//i, showing that intervalley 
scattering in the system is negligible. The resonances 
for fi' < — 3A/2 are due to quasi-bound states in the va- 
lence band of the second constriction [13] . For positive /x' 
(n' ~ 1) we would expect a conductance of order e'^/h, 
because the two constrictions transmit the same valley 
polarization. This is indeed observed for > 1, but re- 
markably enough the conductance remains < 10~^e^/h 
if the arms of the ring support a single propagating mode 
{N = 1, black line in Fig.[2|). We attribute this anomaly 
to current blocking from mismatched valley polarizations 
at the vertices of the hexagonal ring, as explained in Fig. 

m 

In Fig. [4] the conductance G is plotted as a function of 
the flux $ = {S/Sc)^c through the ring, where $c is the 
fiux per unit cell of area S'c = ^V^a"^, and 5' = {So+Si)/2 
is the average of the outer and inner ring areas So and 
Si- The first harmonic frequency of the Fourier spectrum 
shown in Fig. [5] is within a few percent of the expected 
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FIG. 2: Conductance of the hexagonal graphene ring in zero 
magnetic field at fixed fi — 0.05 r ~ A/3 as a function of 
n'/A (with A = 7rr/20 for w = lOVSa). The parameter 
jj,' = — U' is varied by varying U' at fixed fioo- The 
three curves are for different ring widths, corresponding to 
A'^ = 1, 3, 5 propagating modes along the ring. Bottom panel: 
The conductance for = 1 on an expanded (logarithmic) 
scale. 
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FIG. 4: Conductance as a function of magnetic flux through 
the ring. Red solid curves show the case n = n' — 1, fi — fi' = 
0.05 T of identical valley polarizations in both constrictions, 
and blue dashed curves show the case n = —n' — 1, jj. = 
—fi' = 0.05 T of opposite polarizations. The number A*' of 
propagating modes in the ring is varied between the three 
panels. 



FIG. 3: Schematic illustration of the mechanism of current 
blocking by mismatched valley polarizations at a vertex of 
the hexagonal ring with zigzag edges, responsible for the sup- 
pressed conductance when A'^ = 1. Shown is the first Brillouin 
zone at the two sides of the vertex, rotated by 7r/3. Solid ar- 
rows indicate the direction of propagation of the lowest mode 
in one of the two valleys (located near the three red dots in the 
Brillouin zone) . The lowest mode in the other valley (near the 
three black dots) propagates in the opposite direction (dot- 
ted arrow, shown only near one of the black dots for clarity) . 
The current is blocked because the corresponding points in 
the Brillouin zone propagate in opposite directions at the two 
sides of the vertex. 



value <i>Q ^ = e//i, indicating that S accurately represents 
the effective area of the ring. 

The harmonic content of the magnetoconductance 
oscillations is strikingly different when the current is 
blocked {N = 1, top panel in Fig. 2]) and when it is not 
blocked {N = 3,5, lower two panels). On the one hand, 
when the conductance is suppressed below 10""^ e^/ft, the 
magnetoconductance oscillations are nearly sinusoidal, 
almost without higher harmonics. This is as expected 
for transmission through evanescent modes. On the other 
hand, when the conductance is of order e^/h the oscil- 
lations are highly nonsinusoidal, with appreciable higher 
harmonics, as expected for transmission through propa- 
gating modes. 
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FIG. 5: Fourier transform of the magnetoconductance data 
shown in Fig. |4] for A'' = 3 and 5. 



A feature of the nonsinusoidal magnetoconductance os- 
cillations shown in Fig. [H and quantified by the Fourier 
transform in Fig. O is the suppression of the lowest har- 
monic (period A$ = (f>o) in the case of opposite valley 
polarizations in the two constrictions (n ~ —n' — 1). 
This suppression of the fundamental periodicity is dra- 
matic for iV = 3 (top panels in Figs. [J] and [5]), but it 
is also noticeable for N — 5 (lower panels). The second 
(and higher) harmonics, in contrast, are enhanced in the 
case of opposite valley polarizations. 

The suppression of the first harmonic indicates that 
electrons which travel only once along the ring have a 
small probability for intervalley scattering and can there- 
fore not contribute to the conductance when the two 
constrictions have opposite valley polarization. Higher 
harmonics correspond to electrons which travel several 
times along the ring, with a larger probability for in- 
tervalley scattering and therefore a larger probability to 
contribute to the conductance. 



IV. GENERIC GRAPHENE RING 

In order to determine how generic our results for the 
hexagonal ring might be, we have repeated our calcula- 
tions for an approximately circular ring with an irregular 
boundary (see Fig. [S]). We constructed the circular ring 




FIG. 6: Graphene ring with an irregular, approximately circu- 
lar, boundary. Bottom panel: A magnified section of the ring 
with R — 35V3a and W = Wy/Sa used in our simulations, 
which shows the irregularity of the boundary. 



by starting from the hexagonal ring and then keeping 
only lattice sites within an annulus formed by two con- 
centric circles. The electrostatic potential shown in the 
bottom panel of Fig. [1] and the vector potential given by 
Eq. ([3]) remain unchanged, only the length of the cen- 
tral area Lc is replaced now by L'c = ViR^-w^. The 
point contacts still have a zigzag boundary, so the fil- 
tering property at entrance and exit should remain in- 
tact, but the circular ring no longer has zigzag bound- 
aries along its entire perimeter (as it did in the case of 
the hexagonal ring). Results are shown in Figs.[7H9l 

Because the boundaries of the ring are not uniform, 
the number of propagating modes in the ring is not well- 
defined. We still label the data by the same number N 
as in the hexagonal ring, but now this number refers only 
to the number of modes which the ring would support if 
the boundaries were of zigzag form. This uncertainty in 
the definition of N might explain why we now see the 
current blocking not only for = 1 but also for iV = 3 
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FIG. 7: Conductance of the circular graphene ring in zero 
magnetic field at fixed fx — 0.05 r as a function of fi' / A. The 
lower panel shows the data for TV — 1, 3 on an expanded 
(logarithmic) scale. 



(Fig. [71 lower panel) . Only for = 5 do we obtain an 
appreciable current through the ring. The conductance 
for = 5 is still well below the optimal value of 2e^//i, 
but more than 100 times larger than for N = 1,3. 

The harmonic content of the magnetoconductance os- 
cillations (Figs, m [9]) is qualitatively similar as in the 
case of the hexagonal ring. When the current is blocked 
{N — 1,3) we see nearly sinusoidal oscillations, almost 
without higher harmonics. When the current is not 
blocked (N = 5) higher harmonics do appear and the 
lowest harmonic is strongly suppressed when the polarity 
of the valley filters in the two point contacts is opposite. 



V. CONCLUSIONS 

In conclusion, we have identified signatures of valley 
polarization in the magnetoconductance of an Aharonov- 
Bohm interferometer in graphene. The suppression of 
the lowest harmonic of the conductance oscillations that 
appears when the two point contacts have opposite val- 
ley polarity indicates that electrons which have travelled 
only once along the ring preserve their valley polariza- 
tion, and therefore cannot contribute to the conductance. 
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FIG. 8: Conductance of the circular ring as a function of 
magnetic flux. 
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FIG. 9; Fourier transform of the magnetoconductance data 
of Fig. E for iV = 5. 
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The special case of a single-mode hexagonal ring 
shows nearly sinusoidal magnetoconductance oscillations 
around a greatly suppressed average conductance, at- 
tributed to a current blocking effect at a vertex where 
two zigzag edges meet at an angle of 7r/3. This current 
blocking is not a special feature of the hexagonal geome- 
try. Instead, we have found the same current blocking in 
a circular ring with irregularly shaped boundaries. This 
finding is consistent with a recent theory for the bound- 
ary condition of a terminated honeycomb lattice at an ar- 
bitrary crystallographic orientation [ll| . It is found that 
the zigzag boundary condition applies generically for any 
angle (t) ^ Q (mod 7r/3) of the boundary (where (j) — Q 
labels the armchair orientation). The current blocking 



mechanism illustrated in Fig. [3] for the vertex between 
two zigzag boundaries would then apply generically to 
any boundary near an orientation of = (mod tt/B). 
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